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Abstract 



We present calculations of several pp scattering cross sections with potential 
applications at the LHC. Significantly large rates for momentum fraction, 
x, as low as 10~ 7 are obtained, allowing for possible extraction of quark and 
\ gluon densities in the proton and nuclei down to these small x values provided 

ON ■ a detector with good acceptance at maximal rapidities is used. 
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I. INTRODUCTION 



In this paper, we study the measurement of quark and gluon distribution functions inside 
the proton at very small momentum fractions x. We consider several processes in pp collisions 
at a center-of-mass energy of 14 TeV and show that the range of useful measurements 
extends down to about x ~ 10~ 6 for most processes, and even down to x ~ 10" 7 for the 
Drell-Yan process. [] Such measurements require a detector which has sufficient acceptance at 
maximal rapidities. Plans for building a detector of this type (FELIX) are under discussion 
at the LHC []. We show how the parton densities may be determined. The event rates are 
high; to estimate them, we use the CTEQ3M distributions [] to provide an extrapolation 
from the region where measurements currently exist, which is x > 10~ 4 . 

Although our results are obtained for proton-proton interactions, a similar analysis can 
easily be applied to proton-nucleus interactions in the same accelerator. For example, in 
the case of proton-calcium collisions, at say a center-of-mass energy of 63 TeV, the same 
region, down to x ~ 10~ 6 and lower, can also be probed. We show how data from such 
studies would provide information on the parton densities in nuclei at both large and small 
x. Previous data for hard processes on nuclei have been confined to fixed-target energies, 
and so the range of processes for which perturbative QCD (pQCD) calculations are reliable 
is limited. 

At the LHC, the integrated luminosity for a proton-proton luminosity of 10 31 cm -2 s _1 
and a run-time of 10 7 sec/yr is 100 pb~ . However, there is a loss of luminosity in proton- 
nucleus collisions; for protons on calcium, a luminosity of 10 30 cm -2 s _1 with a shorter run- 
ning time is envisaged. This loss of luminosity is to a large extent compensated by increased 
cross sections, which are approximately proportional to the mass number A for large x. At 
small x, the cross sections presumably behave more like A 2 / 3 , but as we will see, the event 
rates are so large that the resulting loss of event-rate relative to proton-proton collisions will 
still enable a lot of physics to be probed. 



II. KINEMATICS AND CROSS SECTIONS 

We consider the following hadronic processes at center of mass energy a/s = 14 TeV with 
s = ( Pl +p 2 ) 2 



p(pl)+p(p 2 ) - 


^jet + 7+ X 


(Jt), 


P(Pl) +P(P2) ' 


^11+ X 


(DY), 


p(pi)+p(p 2 ) - 


-►jet! + jet 2 + X 


(JJ), 


p(pi)+p(p 2 ) ~ 


^Q + Q+ X 


(HQ), 


P(Pi)+P(P2) - 


-> W/Z + X 


(VB), 



where 7 denotes a photon, // a lepton pair, Q(Q) a heavy quark (antiquark) and W, Z are 
the weak vector-bosons. 



1 Similar considerations are valid for the Tevatron collider where due to the smaller ^fs the relevant 
regions in x are scaled up by a factor of ~ 60. 
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To calculate the cross sections, we use the pQCD formalism, with the lowest-order hard 
scattering cross sections found, for example, in Refs. []. For the production of jets, photons 
and heavy quarks, we impose a cut pr,min = 10 GeV, which will keep us in the region where 
perturbative calculations are applicable. Since our main interest is to probe parton densities 
at small x, we will mostly need asymmetric configurations of the momentum fractions, x\ 
and x 2 , of the partons entering the hard scattering. From the kinematic inequality 

. 4 PT,min 
X\X 2 > , 

s 

we deduce that the momentum fractions obey X\X 2 > 2 x 10~ 6 . High-energy data on soft 
hadron production indicate [] that for fixed y ma x — y, where y is the particle rapidity, the 
soft hadronic multiplicity does not increase with energy, although at y — it grows rapidly 
with s. Thus, for the large values of y that we use in this analysis, soft interactions result in 
much smaller underlying event pedestals than for y ~ at LHC energies. Moreover, for 
x substantially larger than x m « n = 2 x 10~ 6 /max(a;i, x 2 ), the counting rates are so high for 
moderate pt of the jets (in processes J7 and J J in expression (|I])) that it would be possible 
to restrict the analysis to the region of sufficiently large x (for one of the partons) so that 
the jets are still produced at large y. Indeed, the expected rates are so high that it would 
be also possible to check the role of the pedestals by taking data at several x bins. Hence, 
our choice of PT,m\n = 10 GeV appears safe. 

Since our interest is only in estimating rates and not in a detailed extraction of parton 
densities from data, it will be sufficient to perform leading-order calculations for most of the 
processes considered. Of course, when actual data become available, it will be necessary to 
fit the parton densities to the data with the aid of theoretical formulae at the best possible 
accuracy, at least next-to-leading order. 

At small x, the gluon density is substantially larger than the quark densities. Since 
at lowest order, the Drell-Yan and vector-boson processes are given by quark- ant iquark 
annihilation, without a gluon-induced subprocess, we will calculate these processes to next- 
to- leading order, with the hard-scattering coefficients in the MS scheme found in Ref . [] . 

For jet, photon and heavy-quark production, we will set the renormalization and factor- 
ization scale /1 to the commonly used value of pt- For Drell-Yan and vector-boson produc- 
tion, we set /i to the pair mass and vector-boson mass, respectively. 



III. JET PLUS PHOTON 

We calculate the cross section for producing a jet and a photon, putting the events in bins 

of Pt and x, where x is the minimum of the momentum fractions, X\ and x 2 , of the incoming 

1 

partons. The x bin is defined by x/A < x < xA, where A = 10 2 ™*" with n^v = 10, while 
the bin size in px is set to 20% of the central value up to px ~ 110 GeV and 40% beyond 
this pt value. These same bin sizes are also used in the other cross sections computed in 
this paper. 

The cross sections shown in Fig. |T| correspond to the kinematic region 

(xi<x 2 <0.8) U (x 2 <xi <0.8). (2) 
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FIG. 1. The cross section ct(x,pt) for Jet + 7 production as a function of x = min(xi, £2) and 
Pt- It is integrated over a pt bin which is 20-40% of the central value (see text), max(xi,X2) < 0.8 
and x/A < min(xi, X2) < xA, where A = 10 1//20 . This choice of A corresponds to 10 bins per decade 
in x. 

One sees that large cross sections (above 10 pb) are obtained for px up to about 100 GeV. 
Combined with a luminosity of 100pb _1 /yr, this gives at least hundreds of events in every 
bin, which will give good statistical precision. The strong fall-off of the curves at their left 
end is a result of approaching the kinematic limit: x > 4p^ min /s. To get a useful number of 
events, it is sufficient for x to be about twice its limit. 

The cross section is dominated by gluon-quark scattering. In Fig. 0, the cross section 
integrated over p T > 10 GeV is split into gq and qq components, and we see that the gq term 
is about an order of magnitude larger over the whole range shown. As before, the fall-off at 
the left is a consequence of the chosen minimum px- Since the incoming quark is typically 
at large x, where its distribution is already known fairly accurately [], photon-jet production 
provides a direct measurement of the gluon density for small x in the range x>2.5xl0~ 6 . 

To illustrate the accuracy that can be achieved in a determination of the gluon density, 
we show in Fig. |^ the CTEQ3M gluon momentum density xG(x, Q 2 = pj) together with the 
statistical errors on measurements that correspond to the cross sections in Fig. |I[ Notice 
that there is sufficient precision not merely to measure the gluon density but also to test its 
evolution. The cross sections presented here are, of course, the result of an extrapolation of 
parton densities from a region where they have been measured to much smaller x values. 
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FIG. 2. The cross section a{x) for Jet + 7 production as a function of x = min(xi,X2)- It is 
integrated over > 10 GeV and the same region in x±,X2 as in the caption for Fig. 111. 



Q= 1 OGeV 




FIG. 3. The gluon momentum distribution with error bars calculated from the Jet + 7 cross section 
and a data sample corresponding to an integrated luminosity of 100 pb _1 . 
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FIG. 4. The cross section a(x, Q) for Drell-Yan as a function of x = min(xi,X2) and Q. It is 
integrated over a Q bin which is 20-40% of the central value and the same region in xi,X2 as in the 
caption for Fig. ffl. 

IV. LEPTON PAIRS 

We next consider lepton pair production from hadrons (the Drell-Yan process). Using 
the same kinematic region defined by expression (Q), we obtain the next-to- leading order 
cross sections shown in Fig. Again, the relatively large cross sections lead to a large 
number of events expected at the LHC. More significantly, by measuring cross sections at 
smaller values of the pair mass Q, we can obtain parton densities at lower values of x than 
in the other processes considered in this paper. The advantage of the Drell-Yan process is 
that one can go to fairly small Q values while still trusting the pQCD formalism. 

The dominant contribution for this process comes from the uu channel; the qg contribu- 
tions are about 30 times smaller than those from qq. Using quark distributions q(x2) which 
are well-determined for large X2, one may then extract antiquark densities q(x) at small x. 
We show in Fig. |5] a plot of the extrapolated CTEQ3M antiquark density xu(x, Q 2 ) with 
the statistical errors based on the cross sections shown in Fig. [|. Note that we consider the 
Drell-Yan process down to relatively small pair masses ~ 5 GeV. This assumes that, in the 
experiment, the detector will be able to suppress the background due to heavy flavor decays. 
Such suppression might be achieved using information from the forward calorimeters as well 
as from microvertex detectors. 
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FIG. 5. The u antiquark momentum distribution with error bars calculated from the Drell-Yan cross 
section and a data sample corresponding to an integrated luminosity of 100 pb _1 . 



V. TWO JETS 

We next consider the cross section for the production of dijets. Fig. |G] shows the cross 
section in bins of px and x, in the region defined by expression (|2]). The cross sections 
obtained are about 100 times larger than all the others computed in this paper. 

In Fig. |7|, we show the x distribution of the cross section. Aside from the total cross 
section (solid curve), we also exhibit the contributions from the different partonic channels: 
gg (dashed), gq (upper dotted), qq (lower dotted) and qq (dot-dashed curve). The largest 
contributions to the cross section come from the gg and gq channels. This process thus 
provides an independent consistency check of the parton distributions, primarily the gluon 
density, obtained from other processes. 



VI. W PRODUCTION 



Production of VF-bosons in the forward direction provides a complementary way to 
measure sea quark distributions at small x. At the same time, comparison of the rates 
of production of W + and W~ bosons would provide a new method to measure the ratio 
r(x, Q 2 ) = ^'q2) at large x [}. (Such a measurement would rely on the presumed equal- 
ity of u and d distributions at small x.) Current determinations olu/d rely on data from 
deuterium targets, where nuclear effects can be large for x > 0.5. Notice that most of the 
current global fits, eg. [], assume that r — > when x — > 1, whereas a perturbative QCD 
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FIG. 6. The cross section ct(x,pt) for dijet production as a function of x = min(xi,X2) and pr- 
Same comments as in the caption for Fig. [I] apply. 





FIG. 7. The cross section a{x) for dijet production as a function of x = min(:Ei, X2). Same 
comments as in the caption for Fig. || apply. 
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FIG. 8. The cross section a(x) for W production as a function of x = x\. It is integrated over 
x/A < x < xA, where A = 10 1 / 20 . This choice of A corresponds to 10 bins per decade in x. Note 
that for this process, x^ = M^/{xs). 

analysis of the leading diagrams for the x — > 1 limit [] suggests that r — > 0.2 in this limit. 
A direct measurement independent of any data requiring knowledge of nuclear effects would 
be valuable. 

At next-to-leading order, we calculate the cross sections for W + and W~ production in 
bins of x. One observes from Fig. [S] that significant rates can be measured up to x ~ 0.8. 
This implies that quark distributions at Q 2 ~ 10 4 GeV 2 can be measured from the sum of 
the cross sections for production of W + and VF~-bosons down to x ~ 5 x 10~ 5 . At the 
same time, one would be able to distinguish between different scenarios for the asymptotic 
behaviour of the u/d ratio in the x — > 1 limit. 



VII. HEAVY QUARKS 

We show in Figs. |^ and |1(] the (x, pt) and x distributions, respectively, for charm quark 
production in the region defined by expression (0). The same general features as in the plots 
for jet plus photon and dijets are observed. 

In Fig. [10], the contributions from the two active partonic channels are also shown: gg 
(dashed curve) and qq (dot-dashed curve). The total cross section is of order 1 to 10 2 pb 
when 2.5 x 10~ 6 < x < 4 x 10 -6 . In this region, the qq contribution dominates, being 
about 10 times larger than that of the gg channel. In the region 10~ 4 < x < 10 -3 , the total 
cross section rises, is of order 10 5 pb and mostly consists of gg channel contributions; the qq 
contribution drops to about 10% of the total. 
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FIG. 9. The cross section a(x,px) for charm quark production as a function of x = min(xi, X2) 
and pr. Same comments as in the caption for Fig. [I] apply. 
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FIG. 10. The cross section a(x) for charm quark production as a function of x = min(:ci, a^)- 
Same comments as in the caption for Fig. ^| apply. 
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Thus, it is possible to extract the antiquark (or quark) density at x ~ 3— 4x 10~ 6 from the 
charm quark cross section providing a cross-check of the antiquark densities obtained from 
Drell-Yan measurements. In this case, the background consisting of gg channel contributions 
is well determined with the use of the gluon density obtained from the jet plus photon or 
dijet cross section as described in Sections [TIT] and |V|. Similarly, it is also possible to extract 
the gluon density at x > 10 -4 from this type of cross section, providing a cross-check of 
gluon densities derived from other measurements. 



VIII. CONCLUSIONS 

We have shown, in several independent ways, how parton distributions in nucleons and 
nuclei (nuclear shadowing) down to x ~ 10~ 7 may be measured from hadron-hadron inter- 
actions at LHC energies. In particular, we show in Figs. |lj] and |12| the regions in energy 
E = (px, Q or M\y) and momentum fraction x where measurements can be made for ex- 
tracting antiquark and gluon densities, respectively, at better than 20% accuracy from the 
different hadronic processes considered, assuming an integrated luminosity of 100 pb~\ 



In Fig. [IT], the region marked with cross-hatch shading shows where the Drell-Yan process 
determines the antiquark densities with better than 20% accuracy. The dashed horizontal 
line is where W production determines these densities. The small triangle with horizontal 
shading shows where charm production determines the antiquark densities. The dotted 
curve shows the minimum value of x as a function of the pair mass Q in the Drell-Yan 
process. 

Similarly, we show in Fig. |T2] the regions where data on the different processes can 
determine the gluon density. The regions marked with vertical, horizontal and slanted lines 
correspond to jet plus photon, charm quark and dijet production, respectively. The dot- 
dashed curve represents the minimum allowed value of a: as a function of the px of the jet 
or heavy quark. 

In our calculations, we have assumed that the usual DGLAP evolution equations apply 
to the parton densities in the very small x region being considered and that CTEQ3M parton 
distributions can be extrapolated outside the region in which they are valid. However, it 
may well be that other physics effects, notably gluon recombination [], are important at 
these small x values. As we have seen, by measuring parton densities over a wide range of 
scales and in several processes, one can, in fact, test whether DGLAP evolution and the 
rest of the factorization method is valid — see, for example, Figs. |3] and |5|. Of course, if the 
measurements fail consistency checks, for example if DGLAP evolution is found to fail, then 
one can no longer say that it is the parton densities that are being measured, but rather 
that one is probing new effects at small x. 
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FIG. 11. A contour plot showing the region in (E, x) where measurements of the antiquark densities 
in the proton can be made to within 20% accuracy, assuming an integrated luminosity of 100 pb _1 . The 
region marked with cross-hatch shading corresponds to the Drell-Yan process. The dashed horizontal 
line is where W production dominates. The dotted curve shows the minimum value of x as a function of 
the pair mass in the Drell-Yan process. The small triangle with horizontal shading shows where charm 
production is dominated by the qq — ► cc subprocess. 
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FIG. 12. A contour plot showing the region (E,x) where measurements of the gluon density in 
the proton can be made to within 20% accuracy, assuming an integrated luminosity of 100 pb _1 . The 
dot-dashed curve represents the minimum allowed value of x as a function of the px of the jet or heavy 
quark. The regions marked with vertical, horizontal and slanted lines correspond to jet plus 7, charm 
quark and dijet production, respectively. 
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